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Applicant's or agent's file reference 
C PH/vdyck/2P 


IMPORTANT NOTIFICATION 


International application No. 
PCT/BE99/00076 


International filing date (day/month/year) 
14 June 1999 (14.06.99) 



1. The following indications appeared on record concerning: 
| | the applicant | [ the inventor X the agent | | the common representative 


Name and Address 

LAND, Addick Adrianus Gosling 

Arnold & Siedsma 

39 Avenue de la Faisanderie 

B-1150 Brussel 

Belgium 


State of Nationality 


State of Residence 


Telephone No. 
32 2 779 44 82 


Facsimile No. 

32 2 779 44 84 


Teleprinter No. 


2. The International Bureau hereby notifies the applicant that the following 
| | the person [ [ the name X the address [ 


change has been recorded concerning: 
| the nationality | | the residence 


Name and Address 

LAND, Addick Adrianus Gosling 
Arnold & Siedsma 
Sweelinckplein 1 
NL-2517 GK The Hague 
Netherlands 


State of Nationality 


State of Residence 


Telephone No. 

31 70 365 4833 


Facsimile No. 

31 70 345 2140 


Teleprinter No. 


3. Further observations, if necessary: 


4. A copy of this notification has been sent to: 
| X| the receiving Office | | the designated Offices concerned 
| | the International Searching Authority the elected Offices concerned 
| X[ the International Preliminary Examining Authority | | other: 





Authorized officer 


The International Bureau of WIPO 


34, chemin des Colombettes 


Sean Taylor 


121 1 Geneva 20, Switzerland 


Facsimile No.: (41-22) 740.14.35 


Telephone No.: (41-22) 338.83.38 
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Assistant Commissioner for Patents 
United States Patent and Trademark 
Office 
Box PCT 

Washington, D.C. 20231 
ETATS-UNIS D'AMERIQUE 

in its capacity as elected Office 


Date of mailing (day/month/year) 
24 August 2000 (24X)8.00) 




International application No. 
PCT/BE99/00076 


Applicant's or agenfs file reference 
C PH/vdyck/2P 


International filing date (day/month/year) 
14 June 1999(14.06.99) 


Priority date (day/month/year) 
16 June 1998 (16.06.98) 


Applicant 

VAN DYCK, Dirk, Ernst, Maria et al 



1. The designated Office is hereby notified of its election made: 

| X | in the demand filed with the International Preliminary Examining Authority on: 

11 January 2000 (11. 01 00) 

| | in a notice effecting later election filed with the International Bureau on: 



2. The election 




not 



made before the expiration of 19 months from the priority date or, where Rule 32 applies, within the time limit under 
Rule 32.2(b). 





Authorized officer 


The International Bureau of WIPO 


34, chemin des Colombettes 


Olivia RANAIVOJAONA 


1211 Geneva 20, Switzerland 




Facsimile No.: (41-22) 740.14.35 


Telephone No.: (41-22) 338.83.38 
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INTERNATIONAL PRELIMINARY EXAMINATION REPORT 

(PCT Article 36 and Rule 70) 



Applicant's or agent's file reference 
C PH/VM81/2p 


_ _ _ See Notification of Transmittal of International 
FOR FURTHER ACTION Preliminary Examination Report (Form PCT/IPEA/416) 


International application No. 
PCT/BE99/00076 


International filing date (day/month/year) 
14/06/1999 


Priority date (day/month/year) 
16/06/1 998 


International Patent Classification (IPC) or national classification and IPC 
H01J37/302 


Applicant 






VAN DYCK, Dirk, Ernst, Maria et al. 







1. This international preliminary examination report has been prepared by this International Preliminary Examining Authority 
and is transmitted to the applicant according to Article 36. 

2. This REPORT consists of a total of 5 sheets, including this cover sheet. 

S This report is also accompanied by ANNEXES, i.e. sheets of the description, claims and/or drawings which have 
been amended and are the basis for this report and/or sheets containing rectifications made before this Authority 
(see Rule 70.16 and Section 607 of the Administrative Instructions under the PCT). 

These annexes consist of a total of 4 sheets. 



3. This report contains indications relating to the following items: 



Basis of the report 



Reasoned statement under Article 35(2) with regard to novelty, inventive step or industrial applicability; 
citations and explanations suporting such statement 

Certain documents cited 

Certain defects in the international application 

Certain observations on the international application 



II 


□ 


III 


□ 


IV 


□ 


V 




VI 


□ 


VII 




VIII 





Date of submission of the demand 
11/01/2000 


Date of completion of this report 
24.11.2000 


Name and mailing address of the international 
preliminary examining authority: 

^v" European Patent Office 
Wh) D-80298 Munich 
^ Tel. +49 89 2399 - 0 Tx: 523656 eornu d 
Fax: +49 89 2399 - 4465 


Authorized officer ^-zr^ 

Van den Berg, G (| |) 

\3 — — 

Telephone No. +49 89 2399 2499 
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INTERNATIONAL PRELIMINARY 
EXAMINATION REPORT 



International application No. PCT/BE99/00076 



I. Basis of the report 

1 . This report has been drawn on the basis of (substitute sheets which have been furnished to the receiving Office in 
response to an invitation under Articfe 14 are referred to in this report as "originally filed" and are not annexed to 
the report since they do not contain amendments (Rules 70. 16 and 70.17).): 
Description, pages: 

1-12 as originally filed 

Claims, No.: 

1-18 as received on 30/1 0/2000 with letter of 30/1 0/2000 

Drawings, sheets: 

1/4-4/4 as originally filed 



2. With regard to the language, all the elements marked above were available or furnished to this Authority in the 
language in which the international application was filed, unless otherwise indicated under this item. 

These elements were available or furnished to this Authority in the following language: , which is: 

□ the language of a translation furnished for the purposes of the international search (under Rule 23.1 (b)). 

□ the language of publication of the international application (under Rule 48.3(b)). 

□ the language of a translation furnished for the purposes of international preliminary examination (under Rule 
55.2 and/or 55.3). 

3. With regard to any nucleotide and/or amino acid sequence disclosed in the international application, the 
international preliminary examination was carried out on the basis of the sequence listing: 

□ contained in the international application in written form. 

□ filed together with the international application in computer readable form. 

□ furnished subsequently to this Authority in written form. 

□ furnished subsequently to this Authority in computer readable form. 

□ The statement that the subsequently furnished written sequence listing does not go beyond the disclosure in 
the international application as filed has been furnished. 

□ The statement that the information recorded in computer readable form is identical to the written sequence 
listing has been furnished. 

4. The amendments have resulted in the cancellation of: 

□ the description, pages: 

□ the claims, Nos.: 
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INTERNATIONAL PRELIMINARY 

EXAMINATION REPORT International application No. PCT/BE99/00076 



□ the drawings, sheets: 

5. □ This report has been established as if (some of) the amendments had not been made, since they have been 
considered to go beyond the disclosure as filed (Rule 70.2(c)): 

(Any replacement sheet containing such amendments must be referred to under item 1 and annexed to this 
report.) 



6. Additional observations, if necessary: 



V. Reasoned statement under Article 35(2) with regard to novelty, inventive step or industrial applicability; 
citations and explanations supporting such statement 

1. Statement 



Novelty (N) 


Yes: 


Claims 


1 - 18 




No: 


Claims 


none 


Inventive step (IS) 


Yes: 


Claims 


1 - 18 




No: 


Claims 


none 


Industrial applicability (IA) 


Yes: 


Claims 


1 - 18 




No: 


Claims 


none 



2. Citations and explanations 
see separate sheet 

VII. Certain defects in the international application 

The following defects in the form or contents of the international application have been noted: 
see separate sheet 

VIII. Certain observations on the international application 

The following observations on the clarity of the claims, description, and drawings or on the question whether the 
claims are fully supported by the description, are made: 
see separate sheet 



: orm PC7/IPEA/^09 fBov.-c !.\/in ?h&ai f !m!v ip°F^ 




INTERNATIONAL PRELIMINARY International application No. PCT/BE99/00076 
EXAMINATION REPORT - SEPARATE SHEET 

To point V : 

Micro- and Nano -Engineering 97. MNE International Conference on Micro- and 
Nanofabrication", Athens, Greece, September 1997, published in Microelectronic 
Engineering, vol. 41-42, pages 195 - 198, March 1998 relates to proximity effects 
correction for electron beam lithography processes by means of an artificial neural 
network. The presence of a coating on a substrate is considered implicit (cf. 
"Introduction"). The document further teaches the training of a neural network by 
inputting a training pattern and an associated exposure, establishing initial weights 
of the neural network, and determining with the neural network the exposure 
required to obtain a desired pattern (cf. page 197: "Proximity neurocorrector"; 
Figure 4). 

1 . (Novelty) Neither the disclosure of this document nor that of the other (even 
more remote) documents cited in the international search report is prejudicial to 
the novelty of the subject-matter of claims 1 - 18. In particular, no method is 
disclosed which allows for predicting the dose of each pattern point of a desired 
pattern using a regularisation parameter and a smearing function (claim 1). The 
subject-matter of claims 1 - 18 would therefore meet the requirement of Article 
33(2) PCT. 

2. (Inventive step) The subject-matter of claims 1-18 does not seem to be 
suggested by the disclosure of any of the documents cited in the international 
search report nor by a combination thereof. There is no incentive for steps of 
determining the smearing function and the precompensated pattern in the way 
claimed (claim 1 ). The subject-matter of claims 1-18 would therefore meet the 
requirement of Article 33(3) PCT. 

3. (Industrial applicability) The subject-matter of claims 1-18 meets the 
requirement of Article 33(4) PCT. 



To point VII : 

A document reflecting the prior art described on pages 1 and 7 is not identified in 
the description (Rule 5.1(a)(ii) PCT). 
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EXAMINATION REPORT - SEPARATE SHEET 



To point VIII : 

The present application does not meet the requirement of Article 6 PCT: 

Claim 17 is not clear in view of the back references to claims 1 - 6 which do not 
relate to the inputting of any pattern into a neural network. Furthermore, claim 17 
mentions a neural network with weighting factors determined as claimed in any of 
the foregoing claims. It does however not become apparent what weighting ~ 
factors are meant because weighting factors are not identified in the foregoing 
claims. 

The features of the claims are not provided with reference signs placed in 
parentheses (Rule 6.2(b) PCT). 
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INTERNATIONAL SEARCH REPORT 

(PCT Article 18 and Rules 43 and 44) 



Applicant's or agent's file reference 

C PH/vdyck/2P 


FOR FURTHER see Notification of Transmittal of International Search Report 

(Form PCT/ISA/220) as well as, where applicable, item 5 below. 

ACTION 


International application No. 

PCT/BE 99/00076 


International filing date (day/month/year) 

14/06/1999 


(Earliest) Priority Date (day/month/year) 

16/06/1998 


Applicant 

VAN DYCK, Dirk, Ernst, Maria et al . 



This International Search Report has been prepared by this International Searching Authority and is transmitted to the applicant 
according to Article 18. A copy is being transmitted to the International Bureau. 

This International Search Report consists of a total of 3 sheets. 

[X] It is also accompanied by a copy of each prior art document cited in this report. 



1 . Basis of the report 

a. With regard to the language, the international search was carried out on the basis of the international application in the 
language in which it was filed, unless otherwise indicated under this item. 

I I the international search was carried out on the basis of a translation of the international application furnished to this 
Authority (Rule 23.1 (b)). 

With regard to any nucleotide and/or amino acid sequence disclosed in the international application, the international search 
was carried out on the basis of the sequence listing : 
| | contained in the international application in written form. 

filed together with the international application in computer readable form. 



b. 



2. 
3. 



□ 
□ 
□ 
□ 

□ 

□ 
□ 



furnished subsequently to this Authority in written form, 
furnished subsequently to this Authority in computer readble form. 

the statement that the subsequently furnished written sequence listing does not go beyond the disclosure in the 
international application as filed has been furnished. 

the statement that the information recorded in computer readable form is identical to the written sequence listing has been 
furnished 

Certain claims were found unsearchable (See Box I). 
Unity of invention is lacking (see Box II). 



4. With regard to the title, 

[~X~| the text is approved as submitted by the applicant. 

| | the text has been established by this Authority to read as follows: 



5. With regard to the abstract, 

|X| the text is approved as submitted by the applicant. 

□ the text has been established, according to Rule 38.2(b), by this Authority as it appears in Box III. The applicant may, 
within one month from the date of mailing of this international search report, submit comments to this Authority. 

6. The figure of the drawings to be published with the abstract is Figure No. 

| | as suggested by the applicant. [X] None of the figures. 

| | because the applicant failed to suggest a figure. 

| | because this figure better characterizes the invention. 
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INTERNATIONAL SEARCH REPORT " ^fc; 



national Application No 

PCT/BE 99/00076 



A. CLASSIFICATION OF SUBJECT MATTER 

IPC 6 H01J37/302 



According to International Patent Classification (IPC) or to both national classification and IPC 



B. FIELDS SEARCHED 



Minimum documentation searched (classification system followed by classification symbols) 

IPC 6 H01J 



Documentation searched other than minimum documentation to the extent that such documents are included in the fields searched 



Electronic data base consulted during the international search (name of data base and, where practical, search terms used) 



C. DOCUMENTS CONSIDERED TO BE RELEVANT 



Category ° 


Citation of document, with indication, where appropriate, of the relevant passages 


Relevant to claim No. 


A 


JEDRASIK P ET AL: "Optimal filtering 


1,2,7,19 




versus regul ari zation methods in the 






Fourier precompensation based proximity 






neurocorrection in electron beam 






1 ithography" 






MICRO- AND NAN0- ENGINEERING 97. MNE 






INTERNATIONAL CONFERENCE ON MICRO- AND 






NANOFABRICATION, ATHENS, GREECE, 15-18 






SEPT. 1997, 






vol. 41-42, pages 195-198, XP004111700 






ISSN 0167-9317, Microelectronic 






Engineering, March 1998, Elsevier, 






Netherl ands 




X 


page 195, right-hand column, line 1 -page 


18 




196, right-hand column, line 4 






page 197, left-hand column, paragraph 2 






-page 198, right-hand column, last line 






-/- 





Further documents are listed in the continuation of box C. 



Patent family members are listed in annex. 



° Special categories of cited documents : 

"A" document defining the general state of the art which is not 

considered to be of particular relevance 
"E" earlier document but published on or after the international 

filing date 

"L" document which may throw doubts on priority clatm(s) or 
which is cited to establish the publication date of another 
citation or other special reason (as specified) 

"O" document referring to an oral disclosure, use, exhibition or 
other means 

"P" document published prior to the international filing date but 
later than the priority date claimed 



T" later document published after the international filing date 
or priority date and not in conflict with the application but 
cited to understand the principle or theory underlying the 
invention 

"X" document of particular relevance; the claimed invention 
cannot be considered novel or cannot be considered to 
involve an inventive step when the document is taken alone 

"Y" document of particular relevance; the claimed invention 

cannot be considered to involve an inventive step when the 
document is combined with one or more other such docu- 
ments, such combination being obvious to a person skilled 
in the art. 

"&" document member of the same patent family 



Date of the actual completion of the international search 

7 October 1999 


Date of mailing of the international search report 

14/10/1999 


Name and mailing address of the ISA 

European Patent Office, P.B. 5818 Patentlaan2 
NL - 2280 HV Rijswijk 
Tel. (+31-70) 340-2040, Tx. 31 651 epo nl, 
Fax: (+31-70) 340-3016 


Authorized officer 

Schaub, G 
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C.(Contlnuatlon) DOCUMENTS CONSIDERED TO BE RELEVANT 



Category * 



Citation of document, with indication, where appropriate, of the relevant passages 



Relevant to claim No. 



JEDRASIK P: "Neural networks application 
for fast, direct correction kernel 
generation for proximity effects 
correction in electron beam lithography" 
MICRO- AND NANOENGINEERING 94. 
INTERNATIONAL CONFERENCE ON MICRO- AND 
NANOFABRICATION, DAVOS, SWITZERLAND, 26-29 
SEPT. 1994, 

vol. 27, no. 1-4, pages 191-194, 
XP004025063 

ISSN 0167-9317, Microelectronic 
Engineering, Feb. 1995, Netherlands 
the whole document 

EP 0 443 249 A (AMERICAN TELEPHONE & ' 
TELEGRAPH) 28 August 1991 (1991-08-28) 
the whole document 

FRYE R C: "ADAPTIVE NEURAL NETWORK 
ALGORITHMS FOR COMPUTING PROXIMITY EFFECT 
CORRECTIONS" 

JOURNAL OF VACUUM SCIENCE AND TECHNOLOGY: 
PART B, 

vol. 9, no. 6, 

1 November 1991 (1991-11-01), pages 
3054-3058, XP000268519 
the whole document 



1,18 



1,18 



1,18 
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Ration on patent family members 
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(national Application No 

CT/BE 99/00076 



Patent document 
cited in search report 



Publication 
date 



Patent family 
member(s) 



Publication 
date 



EP 0443249 



28-08-1991 



DE 
DE 
JP 
US 



69030078 D 
69030078 T 
6097057 A 
5283746 A 



10-04-1997 

07- 08-1997 

08- 04-1994 
01-02-1994 
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REQUEST 



The undersigned requests that the present 

international application be processed 
according to the Patent Cooperation Treaty. 



:eiving Office use only 

PC T f^E 9 9 / 0 0 0 7 6 

International Application No. 



14 JUIN199 

Internationa] Filing Date 



I (H -06- f99 9 ) 



flO/BE - INTERNATIONAL APPLICATION 



Name of receiving Office and "PCT International Application" 



Applicant's or agent's file reference c PH / vdyck/ 2p 
(if desired) (12 characters maximum) 



Box No. I TITLE OF INVENTION 

IVIethod and device for correcting of the proximityeff ec ts 



Box No. II APPLICANT 



Name and address: (Family name followed by given name; for a legal entity, full official 
designation. The address must include postal code and name of country. The country of the 
address indicated in this Box is the applicant s State (that is, country) of residence if no State 
of resiaence is indicated below.) 

VAN DYCK, Dirk Ernst Maria 
Kleine Grippe 37 
B-2630 AARTSELAAR 



State (th at is, country) of nationality: 



This person is applicant 
for the purposes of: 



This person is also inventor. 



Telephone No. 



Facsimile No. 



Teleprinter No. 



all designated 
States 



□ all designated States except 
the United States of America 



State (that i s^ cou ntry) of residence: 



□ the United States 
of, 



f America only 



□ the States indicated in 
the Supplemental Box 



Box No. IE FURTHER APPLICANT(S) AND/OR (FURTHER) EWENTOR(S) 



Name and address: (Family name followed by given name; for a legal entity, full official 
designation. The address must include postal code and name of country. The country of the 
adaress indicated in this Box is the applicant s State (that is, country) of residence if no State 
of residence is indicated below.) J 

JEDRASIK , Piotr Tomasz 
fcFniversitair Centrum Antwerpen 
Universiteit van Antwerpen 
Departement Natuurkunde 
Groenerborgerlaan 1 71 

R-POPO AMTWKRPFA' 



BELGIUM 



State (that is, country) of nationality: 



This person is applicant 
for the purposes of: 



This person is: 

I I applicant only 

\XX applicant and inventor 



j inventor only (If this check-box 
is marked do not fill in below.) 



State (tha tjs, c ountry) of residence: 



□ all designated I | all designated States except [ ] the United States I — I the States indicated in 

States I 1 the United States of America tod of America only | | the Su^emoS] Box 

□ 

Further applicants and/or (further) inventors are indicated on a continuation sheet- 



Box No. IV AGENT OR COMMON REPRESENTATIVE; OR ADDRESS FOR CORRESPONDENCE 



Tne person identified below is hereby/has been appointed to act on behalf 
of the applicants) before the competent International Authorities as: 



\ck agent | | common representative 



Name and address: (Family name followed by given name; for a legal entity, full official 
designati on. The address must include postal code and name of country.) 



LAND, Addick Adrianus Gosling 

ARNOLD & SIEDSMA 

39, Avenue de la Faisanderie 

B-1150 BRUSSEL 

BELGIUM 



Telephone No. 

+32-2 779 44 82 



Facsimile No. 

+32 2 779 44 84 



Teleprinter No. 



□ Address for correspondence: Mark this check-box where no agent or common representative is/has been appointed and the 
space above is used instead to indicate a special address to which correspondence s hould be sent 
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Sheet No. 



Box No.V 



DESIGNATK 



The following designations ar< 
Regional Patent 

0 



_ OF5 



STATES 



iy made under Rule 4.9(a) (mark the applicable^ck-boxes; at least one must be marked) 



EA 



El EP 



IS 



ARIPO Patent: GH Ghana^ GM Gambia, KE Kenya, LS Lesotho, MW Malawi, SD Sudan, SZ Swaziland, UG Uganda. 
Z\\ Zimbabwe, and any other State which is a Contracting State of the Harare Protocol and of the PCT u sanaa, 
Eurasian Patent: AM Armenia, AZ Azerbaijan, BY Belarus, KG Kyrgyzstan, KZ Kazakhstan, MD Republic of 
Moldova, RU Russian Federation, TJ Tajikistan, TM Turkmenistan, and any other State^chba CMt«2K State 
of the Eurasian Patent Convention and of the PCT s aiaic 

European Patent: AT Austria, BE Belgium, CH and LI Switzerland and Liechtenstein, CY Cyprus, DE Germany 
DK Denmark, ES Spain, FI Finland, FR France, GB United Kingdom, GR Greece, IE Ireland, IT Italy, LU Luxembourg 
MC Monaco, NL Netherlands, PT Portugal, SE Sweden, and any other State which is a Contracting State of the European 
Patent Convention and of the PCT *^ 

° A 2 a V K atC ^5u r F ? urk i2?j5?°.' W i*™* CF Centra] African Republic, CG Congo, CI C6te d'lvoire, CM Cameroon, 
GA Gabon, GN Guinea, GW Guinea-Bissau, ML Mali, MR Mauritania, NE Niger, SN Senegal, TD Chad, TG Togo and 
any other State which js a member State of OAPI and a Contracting State of the PCT Of otter land of protection or *ea\ment 
desired, specify on doaedline) 



National Patent (if other kind of protection or treatment desired, specify on dotted line): 



m 
m 

E 
13 
El 

m 
m 

m 
m 
is 
w 
w 

E3 
E 

m 

HI 

m 

E 
El 
El 
E 
SI 
E 
12 
E 
E] 
IS 
El 
0 
E 
IS 

H 

EJ 

E 
E 



AL Albania 

AM Armenia 

AT Austria 

AU Australia 

AZ Azerbaijan 

BA Bosnia and Herzegovina 

BB Barbados 

BG Bulgaria 

BR Brazil J , 

BY Belarus 

CA Canada 

CH and LI Switzerland and Liechtenstein 

CN China 

Cuba 

Czech Republic 

Germany 

Denmark 

Estonia 

Spain 

Finland 

United Kingdom 
Grenada 

Georgia 

Ghana 

GM Gambia 

HR Croatia 

Hungary * 

Indonesia 

Israel 

India 

Iceland 

Japan 

Kenya 

Kyrgyzstan 

Democratic People's Republic of Korea 




Lesotho 

Lithuania 
Luxembourg 
Latvia 

MD Republic of Moldova 

MG Madagascar 

MK The former Yugoslav Republic of Macedonia 



CU 

cz 

DE 

DK 

EE 

ES 

FI 

GB 

GD 

GE 

GH 



HU 

ID 

EL 

IN 

IS 

JP 

KE 

KG 

KP 



m 

12 
EI 
13 

m 

El 

Q 
EI 

or 
ca 
b 

0 
El 
Si 
□ 



MN Mongolia ' 

MW Malawi 

MX Mexico 

NO Norway 

New Zealand 

Poland 

Portugal , 

Romania 

Russian Federation 

Sudan 
Sweden 
Singapore 

Slovenia 

Slovakia 

Sierra Leone 

Tajikistan 

Turkmenistan 

Turkey 

Trinidad and Tobago . . . 

Ukraine 

Uganda 

United States of America . 



NZ 
PL 
PT 
RO 
RU 
SD 
SE 
SG 
SI 
SK 
SL 
TJ 
TM 
TR 
TT 
UA 
UG 
US 



KR 
KZ 
LC 
LK 
LR 



13 UZ Uzbekistan 
□ VN Viet Nam . 
El YU Yugoslavia 
S ZW Zimbabwe . 

Check-boxes reserved for designating States (for the pu 



<-iicc*-doxcs rcscrvca ior designating Mates (tor the purposes of 

Republic of Korea a national patent) which have become party to the PCT after 

Kazakhstan issuance of this sheet: 

Saint Lucia □ 
Sri Lanka □ 
Liberia □ 



Precautionary Designation Statement: In addition to the designations made above, the applicant also makes under Rule 4 9(b) all other 
designations which would be permitted under the PCT except any designation(s) indicated in the Supplemental Box as being excluded 
from the scope of this statement Hie applicant declares that those additional designations are subject to confirmation and that any 
designation whicn is not confirmed before the expiration of 15 months from the priority date is to be regarded as withdrawn by the applicant I 
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WERKWIJZE EN INRICHTING VOOR HET 
CORRIGEREN VAN NABI JHEIDSEFFECTEN 

De onderhavige uitvinding heeft betrekking op 
een werkwijze en inrichting voor het bepalen van de voor 
verkrijging van een gewenst patroon in een deklaag op een 
ondergrond per patroonpositie benodigde blootstellingsdo- 
5 sis van een elektronenbundel . 

Bij het vervaardigen van de nieuwste generaties 
van geintegreerde schakelingen wordt bij voorkeur gebruik 
gemaakt van gef ocusseerde elektronenbundels in lithogra- 
fische processen in plaats van gebruik te maken van de 

10 gebruikelijke optische lithograf ie- technieken, aangezien 
deze laatste technieken beperkingen met betrekking tot 
resolutie als gevolg van diffractie van het gebruikte 
laserlicht ondervinden. De resolutie van de geintegreerde 
schakeling die verkregen wordt met een dergelijke elek- 

15 tronenbundellithograf ie is groter, doch wordt beperkt 
door verstrooiing van de elektronen in de deklaag. Er 
zijn werkwijzen bekend van het minimaliseren van ver- 
strooiingsef f ecten of van het vooraf hiervoor compenseren 
en het daardoor verbeteren van de resolutie van de ver- 

20 kregen geintegreerde schakelingen. 

De bekende werkwijzen hebben echter als bezwaar 
dat verstrooiingseff ecten zelf slechts in beperkte mate 
te minimaliseren zijn, terwijl een compensatie vooraf 
volgens de bekende werkwijze veel berekeningen vergt en 

25 derhalve een lange berekeningsti j d nodig heeft. Voor het 
vervaardigen van geintegreerde schakelingen dient bij- 
voorbeeld een zeer groot aantal patroonpunten, vaak in de 
orde van groot te van 10 10 patroonpunten, "geschreven" te 
worden, terwijl het aantal hiervoor benodigde bereke- 

3 0 ningen een veelvoud daarvan bedraagt . Hierdoor is een 
nagenoeg onvertraagde ("real time") precompensatie voor 
de versmeringseff ecten niet uitvoerbaar. 
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Het doel van de onderhavige uitvinding is dit 
bezwaar te ondervangen en tevens aanvullende voordelen te 
verschaf f en. 

De onderhavige uitvinding heeft derhalve be- 
5 trekking op een werkwijze voor het bepalen van de voor 
verkrijging van een gewenst patroon in een deklaag op een 
ondergrond per patroonpositie benodigd geprecompenseerd 
patroon van blootstellingsdoses van een elektronenbundel , 
omvattende : 

10 - het bepalen van de versmeringsf unctie van de 

elektronenbundel ; 

- het met de versmeringsf unctie en het gewenst 
patroon bepalen van het geprecompenseerde patroon, 
waarbij het bepalen zodanig wordt uitgevoerd, dat de 

15 blootstellingsdoses vrijwel uitsluitend positieve waarden 
bevatten en dat de blootstellingsdoses ten opzichte van 
elkaar ten minste in enige mate glad zijn. 

Aangezien een negatieve waarde voor de bloot- 
stellingsdosis van een elektronenbundel geen fysische 

20 betekenis heeft en derhalve niet realiseerbaar is, wordt 
de bepaling van de blootstellingsdoses van het geprecom- 
penseerde patroon zodanig uitgevoerd dat deze (vrijwel) 
uitsluitend positieve waarden aannemen. Bovendien wordt 
een gladde oplossing verkregen aangezien hevige oscilla- 

25 ties in de versmeringsf unctie geen fysische grondslag 
hebben, doch slechts veroorzaakt worden. door mathema- 
tische instabiliteit van de berekeningen . 

In een voorkeursuitvoeringsvorm van de uitvin- 
ding omvat de werkwijze de stappen: 

30 a) het schatten van een regularisatie parame- 

ter; 

b) het met alle patroonpunten van het gewenst 
patroon behalve een bepaald patroonpunt bepalen van een 
geprecompenseerd patroon ; 
3 5 c) het met de versmeringsf unctie opnieuw ver- 

smeren van het geprecompenseerd patroon voor het voor- 
spellen van de dosis van het bepaalde patroonpunt; 
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d) het voor elk patroonpunt herhalen van stap- 
pen b en c; 

e) het herhalen van stappen a tot en met d met 
aangepaste regularisatie parameter totdat een eindwaarde 
van een regularisatie parameter bereikt is; 

f) het bepalen van het geprecompenseerde pa- 
troon met de eindwaarde van de regularisatie parameter. 

Volgens een verdere uitvoeringsvorm van de 
uitvinding omvat stap b) de volgende iteratieve bepaling; 

d u) = d^-v + OCK+XBiD) ) -i J Rr*r (I - 1) r a) = a-Kd U) 



met d (0> = 0 en r (0) = a 

waarbij a een vector met als elementen de doses van het 
gewenste patroon, d een vector met de blootstellingsdoses 

15 van het geprecompenseerde patroon, K de versmeringsf unc- 
tie in matrixvorm is, K* de hermitisch geconjugeerde van 
de versmeringsf unctie K is, B een operator, en X een 
regularisat ieparameter . 

Volgens een verdere uitvoeringsvorm van de 

20 uitvinding is operator B als volgt gedef inieerd : 



B(D) (-^-)in(-^-) 



waarin de sommatie plaatsvindt over alle patroonpunten, d. 
het i e . element van de vector d is, en d tot de sommatie over 
25 alle elementen van de vector d voorstelt. 

Volgens een verdere uitvoeringsvorm van de 
uitvinding is in bovengenoemde stap e) de eindwaarde van 
de regularisatie parameter de regularisatie parameter 
waarbij 

N 



30 
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waarin N het totaal aantal patroonpunten, a een vector 
met als elementen de doses van het gewenste patroon, d 
een vector met de blootstellingsdoses van het geprecom- 
penseerde patroon, en K de versmeringsf unctie in matrix- 
5 vorm is. 

Volgens een verdere uitvoeringsvorm van de uit- 
vinding is in stap e) de eindwaarde van de regularisatie 
parameter de regularisatie parameter waarbij 

it£ (a* " l^'wiuttt) 

iV Jt»l 

10 waarin N het totaal aantal patroonpunten, a een vector 
met als elementen de doses van het gewenste patroon, d 
een vector met de blootstellingsdoses van het geprecom- 
penseerde patroon, K de versmeringsf unctie in matrixvorm 
en w^ gedefinieerd is als: 

Wkk a) = [ ]2 

15 met a^ de elementen van de matrix 

A = K(K T K+XL(D) T L(D) )- x K T en L de discrete Laplace -trans - 
formatie, minimaal is. 

Volgens een verdere uitvoeringsvorm van de 
uitvinding worden na stap e) de stap uitgevoerd van het 

2 0 trainen van een neuraal netwerk met behulp van een of 

meer gewenste eerste patronen en de daarbij behorende 
geprecompenseerde patronen. 

Volgens een verdere uitvoeringsvorm van de 
uitvinding is met het getrainde neurale netwerk het 
25 geprecompenseerde patroon behorend bij een tweede gewenst 
patroon te bepalen, waarbij in een verdere uitvoerings- 
vorm het eerste gewenste patroon een relatief eenvoudig 
trainingspatroon is en het tweede gewenste patroon het 
deelpatroon van een geintegreerde schakeling is en waar- 

3 0 bij in een verdere uitvoeringsvorm twee of meer deelpa- 
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tronen tot een samengesteld patroon van de geintegreerde 
schakeling te combineren zi jn . 

Door een bekend eerste gewenst patroon, dat bij 
voorkeur eenvoudig is, het bijbehorende geprecompenseerde 
5 patroon van blootstellingsdosis te bepalen en vervolgens 
het verband tussen de weegfactoren van een neuraal net- 
werk vast te leggen, wordt verzekerd dat voor een tweede 
gewenst patroon, dat gecompliceerder kan zijn, het ver- 
krijgen van het verband tussen dit patroon en de daarbij 

10 behorende blootstellingsdosis op zeer efficiente en 

snelle wijze door het neurale netwerk wordt bepaald. Het 
eerste patroon is een over het algemeen relatief eenvou- 
dig trainingspatroon, terwijl het tweede patroon, bij- 
voorbeeld het patroon van een zeer gecompliceerde , 

15 geintegreerde schakeling is. 

In een voorkeursvorm van de uitvinding is het 
bovengenoemde neurale netwerk in hardware geimplemen- 
teerd, waardoor het bepalen van het verband tussen een 
patroon en de daardoor bijbehorende blootstellingsdosis 

20 op snellere wijze wordt bepaald, bij voorkeur binnen 
6 0 ns voor elk patroonpunt en voor een patroon van 10 10 
binnen de 10 minuten. 

Volgens een voorkeursvorm van de uitvinding is 
de versmeringsfunctie opgebouwd uit ten minste twee 

25 Gaussische functies, waaraan eventueel een exponentiele 
functie is toegevoegd. Parameters van de Gaussische 
functies en eventueel de exponentiele functie zijn te 
bepalen door middel van statistische simulatie van het 
systeem van elektronenbundelzendapparatuur en de betref- 

3 0 fende deklaag en ondergrond van de te vervaardigen 
geintegreerde schakeling. 

In een andere uitvoeringsvorm van de uitvinding 
worden parameters bepaald door metingen te verrichten aan 
het systeem van elektronenbundelzendapparatuur en de 

35 betref fende deklaag met ondergrond. 

De onderhavige uitvinding heeft tevens betrek- 
king op een inrichting voor het bepalen van de voor het 
verkrijgen van een gewenst patroon in een deklaag op een 
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ondergrond per patroonpositie benodigde blootstellingsdo- 
sis van een elektronenbundel , omvattende een elektro- 
nische schakeling voor het implementeren van het bovenge- 
noemde neurale netwerk met weegf actoren die op de boven- 
5 genoemde wijze zijn bepaald. 

De uitvinding zal hierna aan de hand van een 
voorkeursuitvoeringsvorm daarvan verduideli jkt worden, 
waarbij gebruik wordt gemaakt van de bijgevoegde teke- 
ningen , waarin : 

10 - figuur 1 een schematisch overzicht is van een 

voorkeursvorm van een inrichting volgens de uitvinding; 

- figuren 2a-2c een schematische overzicht 
geven van de bepaling van een geprecompenseerde patroon 
van 3x3 patroonpunten; 

15 - figuur 3 een gewenst trainingspatroon van 

256x256 patroonpunten weergeeft; 

- figuur 4 het trainingspatroon van figuur 3 na 
versmering weergeeft; 

- figuur 5 een graf iek weergeeft waarin voor 
20 het trainingspatroon van figuur 3 de voorspellingsf out 

als functie van de gekozen regularisatieparameter is 
uitgezet ; 

- figuur 6 het trainingspatroon van figuur 3 na 
precompensatie weergeeft; 

25 - figuur 7 het geprecompenseerde patroon van 

figuur 6 na versmering weergeeft; en 

- figuur 8 een schematische voorstelling van 
een neuraal netwerk voor de bepaling van geprecompenseer- 
de patronen weergeeft. 

3 0 In een opstelling van apparatuur voor het 

zenden van een elektronenbundel en een te behandelen 
ondergrond 1 met deklaag 2 wordt een bundel elektronen 3 
gericht op een positie of patroonpunt van een deklaag 2 
op een ondergrond 1. De interactie van de invallende 

35 elektronenstraal 3 met de deklaag of resistfilm 2 en de 
onderlaag of het substraat 1 heeft een verstrooiing van 
de elektronen in de deklaag 2 tot gevolg hetgeen versme- 
rings- ofwel nabi jheidsef f ecten veroorzaakt. Wanneer 
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bijvoorbeeld een primair elektron de deklaag binnen- 
dringt, wordt een deel van zijn energie overgebracht op 
elektronen van de atomen van de deklaag, hetgeen ioni- 
satie of excitatie daarvan veroorzaakt . Een botsing 
5 tussen elektronen met een grote overdracht van energie 
genereert een secundair elektron, dat in het algemeen een 
verplaatsingsrichting loodrecht op die van een primair 
elektron heeft. 

Meer in het algemeen hebben versmeringsef f ecten 

10 in elektronenbundellithograf ie betrekking op het proces 
waardoor de resolutie van het blootgestelde patroon 
verlaagd wordt door primaire elektronenverstrooiing 
(forward scattering) en secundaire elektronopwekking 
(backward scattering) in de deklaag en de ondergrond van 

15 een te vervaardigen geintegreerde schakeling. Scherpe 
kenmerken zoals hoeken in het gewenste patroon worden 
afgerond, lijndiktes en tussenruimtes worden aangepast en 
in bepaalde extreme gevallen verdwijnen zelfs sommige 
kenmerken geheel of worden zij op incorrecte wijze samen- 

20 gevoegd met naastliggende kenmerken. 

De versmeringsef f ecten of nabi jheidsef f ecten 
zijn te beschrijven door een versmeringsf unctie , die het 
verband weergeeft tussen enerzijds de blootstellingsdoses 
van een bepaald patroonpunt van een te vervaardigen 

2 5 patroon in de deklaag en anderzijds de door dit patroon- 

punt en door naburige patroonpunten daadwerkeli jk geab- 
sorbeerde doses. Het effect van de versmering wordt der- 
halve vastgelegd in de versmeringsf unctie . 

Ervan uitgaande dat de blootstelling en versme- 

3 0 ring lineair en spatieel invariant zijn en dat voor een 

numerieke oplossing een discrete representatie de voor- 
keur heeft kan het bovenstaande als volgt in matrixvorm 
tot uiting worden gebracht : A = KD, waarin A een kolom- 
vector is waarvan elk element a ± de totale energiedosis is 
35 die daadwerkeli j k in het bijbehorende patroonpunt is 

geabsorbeerd, K een versmeringsmatrix is waarvan elk mn- 
de element het gedeelte van de energiedosis is dat geab- 
sorbeerd is in patroonpunt m vanaf een eenheidsbloot- 
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stellingsdosis die wordt geleverd aan patroonpunt n, en D 
een kolomvector is die is opgebouwd uit elementen d i die 
de blootstellingsdoses representeren die per patroonpunt 
door de elektronenbundelapparatuur wordt afgegeven. 
5 Aangezien het versmeringsef feet onvermi j deli jk is, is het 
zaak om de blootstellingsdoses d i van de verschillende pa- 
troonpunten zodanig aan te passen dat de daadwerkeli jk in 
een patroonpunt geabsorbeerde dosis a i zodanig is dat toch 
het gewenste patroon verkregen wordt . 

10 Deze zogenoemde precompensatie van de bloot- 

stellingsdosis van de elektronenbundel kan volgens de 
huidige stand van de techniek worden uitgevoerd door 
bepaling van de inverse van de versmeringsmatrix K. Er 
zijn vele manieren om een matrix in het algemeen te 

15 inverteren. Deze manieren houden vaak echter geen reke- 
ning met fysieke beperkingen, zoals in dit geval bijvoor- 
beeld die van de elektronenbundel zendapparatuur . Zo zijn 
er bijvoorbeeld geen negatieve blootstellingsdoses moge- 
lijk. Bovendien is een bezwaar van dergelijke inversieme- 

20 thodes dat de geinverteerde matrix veel oscillaties 

vertoont . Bovendien dient voor inversie van de versme- 
ringsmatrix voor een deelpatroon van bijvoorbeeld 256x256 
patroonpunt en al de inversie van een versmeringsmatrix 
met dimensies van 65536x65536 berekend te worden, hetgeen 

25 een enorme hoeveelheid rekentijd vergt . 

In figuren 2a t/m 2c is een gewenst patroon (A) 
weergegeven, Het patroon is opgebouwd uit een 9-tal 
patroonpunten a x waarbij i varieert van 1 tot 9. Dit 
gewenst patroon dient geprecompenseerd te worden teneinde 

3 0 na blootstelling aan de versmerende elektronenbundel het 
gewenste patroon te kunnen verschaf fen, dat wil zeggen 
dat de waarden van d., met i varierend van 1 tot 9, be- 
paald dienen te worden. 

Allereerst wordt het geprecompenseerde patroon 

3 5 bepaald gebruik makend van de doses a ± met i van 2 tot 9, 
waarbij derhalve patroonpunt 1 buiten beschouwing wordt 
gelaten (figuur 2a) . Dit geprecompenseerde patroon wordt 
bepaald aan de hand van de volgende uitdrukking: 
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d d) = d^-v + UCK+XBiD) ) -ijrr (i - 1} r (i) = a-*Td (i) 
met d (0 > = 0 en r< 0) = a 

waarbij a een vector met als elementen de doses van het 
gewenste patroon, d een vector met de blootstellingsdoses 
5 van het geprecompenseerde patroon, K de versmeringsf unc- 
tie in matrixvorm is, K* de hermit isch geconjugeerde van 
de versmeringsf unctie K is, B een operator is, en A een 
regularisatieparameter is. De waarde van de regularisa- 
tieparameter kan willekeurig gekozen worden, in dit geval 
10 bijvoorbeeld k - 0 . 

De operator B legt een beperking op en kan als 
volgt gedefinieerd worden: 

B(D) (J^)ln(-^) 

i "tot "tot 

15 waarin de sommatie plaatsvindt over alle patroonpunten, d. 

het i e element van de vector d is, en d tot de sommatie over 

alle elementen van de vector d voorstelt. 

Het aldus bepaalde geprecompenseerde patroon 

wordt vervolgens aan de hand van de bekende versmerings- 
20 functie opnieuw versmeerd, waarmee de voorspelde dosis Kd 

van patroonpunt 1 wordt bepaald. 

Vervolgens wordt de bovenstaande procedure 

achtereenvolgend herhaald (figuren 2b en 2c) voor het 

tweede tot en met negende patroonpunt (i=2,..,9), waarbij 
25 telkens alle patroonpunten met uitzondering van een 

patroonpunt gebruikt worden. 

Aan de hand van de bovenstaande resultaten 

wordt de kleinste kwadraten voorspellingsf out over alle 

patroonpunten bepaald, hetgeen later nader zal worden 
30 verklaard. 

Vervolgens wordt de bovenstaande procedure 
herhaald met andere waarden voor de regularisatiepara- 
meter A, bijvoorbeeld A 2 = 0,001, X= 0,002, etc. Gekozen 
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wordt uiteindelijk de regularisatieparameter waarbij de 
kleinste kwadraten voorspellingsf out over alle patroon- 
punten minimaal is. Deze regularisatieparameter wordt dan 
gekozen als de optimale regularisatieparameter A . Aan de 
5 hand van deze optimale regularisatieparameter k wordt 

opt 

dan het uiteindeli jke geprecompenseerde patroon bepaald. 

Hiertoe wordt het minimum bepaald van de uit- 

drukking : 
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waarin N het totaal aantal patroonpunten, a een vector 
waarin N het totaal aantal patroonpunten, a een vector 
met als elementen de doses van het gewenste patroon, d 
een vector met de blootstellingsdoses van het geprecom- 
15 penseerde patroon, K de versmeringsf unctie in matrixvorm 
en w^ gedefinieerd is als: 

Wkk a) = [ 1 ~^ U) ]2 

met a^ de elementen van de matrix A = K (K T K+XL (D) T L (D) ) _1 K T 
en L de Laplace-operator. 

De versmeringsf unctie als gevolg van forward 

20 scattering en backward scattering van de elektronen van 
de elektronenbundel kan op verschillende wijzen worden 
bepaald. Zij kan worden bepaald aan de hand van metingen 
van de impulsresponsie van de apparatuur voor het zenden 
van de elektronenbundel op een testobject. De versme- 

25 ringsfunctie kan ook worden bepaald met behulp van diver- 
se. Monte Carlo technieken. In de eerstgenoemde bepalings- 
wijze wordt rekening gehouden met alle fysische aspecten 
van de gebruikte apparatuur. In de later genoemde bepa- 
lingswijze wordt slechts een model van de werkeli jkheid 

30 gebruikt, de bepaling is echter wel eenvoudiger uit te 
voeren zonder uitgebreide metingen te vergen. 
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Als benadering voor de op een van de boven 
beschreven wijzen bepaalde versmeringsf uncties zijn bij 
voorkeur Gaussische f uncties te gebruiken. De versme- 
ringsf unctie wordt in dit geval bijvoorbeeld "gefit" met 
5 een scattering fit model van een dubbele Gaussische 
functie (voor zowel forward als backward scattering 
eigenschappen van de elektronen) , een driedubbele Gaus- 
sische functie of een dubbele Gaussische functie met een 
dalende exponentiele functie. De keuze van het scattering 

10 fit model is afhankelijk van de afmetingen van de te 

onderscheiden onderdelen van het testobject (resolutie) . 
Bij afmetingen kleiner dan 10 0 nm wordt de keuze hiervan 
kritisch; bij dergelijke kleine afmetingen hebben de 
driedubbele Gaussische functies of dubbele Gaussische 

15 functies met dalende exponentiele functie de voorkeur. 

Een versmeringsfunctie met dubbele Gaussische functies is 
met 3 parameters te beschrijven, terwijl de overige twee 
genoemde scattering fit modellen met 4 parameters te 
beschrijven zijn, hetgeen een grote reductie in de hoe- 

20 yeelheid te verwerken gegevens impliceert. 

In figuur 3 is een gewenst patroon van 256x256 
patroonpunten weergegeven. Na versmering met een versme- 
ringsfunctie in de vorm van een dubbele Gaussische func- 
tie met or= 50 nm, £= 3,45 en n = 1,36 levert het versmeer- 

25 de patroon van figuur 4. Duidelijk zichtbaar is, dat veel 
detail in het patroon verloren is gegaan, hetgeen een 
beperking betekent van de te bereiken resolutie van het 
te vervaardigen patroon. Toepassing van de werkwijze 
volgens de uitvinding levert een optimale regularisatie- 

30 parameter van A opt = 0,07042, hetgeen in figuur 5 is aange- 
geven, waarin bij deze waarde van X de fout in het pa- 
troon minimaal is. Het met deze waarde van X berekende 
geprecompenseerde patroon en het bijbehorende versmeerde 
patroon zijn weergegeven in respectieveli jk figuren 6 en 

35 7. Uit vergelijking van de resultaten van figuur 7 met 
die van figuur 3 blijkt, dat het precompensatie van het 
patroon een gewenst patroon een versmeerd patroon ople- 
vert waarvan de resolutie sterk is verbeterd. Hiermee 
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zijn te onderscheiden onderdelen met afmetingen van 
minder dan 100 nm, bijvoorbeeld in een geintegreerde 
schakeling, realiseerbaar . Een vergelijking van de resul- 
taten van de hierin beschreven werkwijze met die van 
5 andere correctiemethoden is opgenomen in tabel 1. De 
foutmaat van de correctiemethoden is hier gedefinieerd 
als de sommatie van het verschil tussen de berekende 
blootstellingsdoses en de ideale geprecompenseerde bloot- 
stellingsdoses gedeeld door het aantal patroonpunten. 

10 



Correct iemethode 


foutmaat in % 


Ongecorrigeerd 


10,2 % 


Af kappen 


10,2 % 


Verschuiven en opschalen 


12,2 % 


Onderhavige methode 


4,9 % 



Uit het bovenstaande blijkt, dat de onderhavige methode 
van bepaling van een geprecompenseerd patroon in verge- 
lijking tot de overige gangbare methoden veruit de klein- 
2 0 ste foutmaat oplevert. 

Het geprecompenseerde patroon en het gewenste 
patroon worden vervolgens gebruikt als trainingsset of 
trainingspatronen voor een neuraal netwerk. Een deel van 
een dergelijk netwerk is schematisch weergegeven in 

2 5 figuur 8 en wordt gerepresenteerd door de uitdrukking 

9 

a i = E w ij h ±j (x) 

dat wil zeggen dat dosis a. uitgedrukt is in een verzame- 
ling van 9 basisf uncties h... , in dit geval radiele func- 
ties . 

3 0 Na training van het neurale netwerk kan voor 

een ander willekeurig gewenst patroon op zeer snelle 
wijze een geprecompenseerd patroon worden bepaald. Een 
willekeurig patroon kan bijvoorbeeld een patroon van 512 
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bij 512 patroonpunten zijn die een deelpatroon van een 
gelntegreerde schakeling vormen. Verscheidene deelpatro- 
nen worden vervolgens gecombineerd (geclusterd) tot een 
patroon dat de gehele gelntegreerde schakeling of ten 
5 minste een deel daarvan omvat . 

Het boven beschreven neurale netwerk kan ge- 
implementeerd worden in hardware, en bij voorkeur in 
analoge hardware aangezien de rekensnelheid van op derge- 
lijke wijze gelmplementeerde neurale netwerken zeer groot 

10 is. Zo bedraagt de rekentijd voor precompensatie van een 
patroon minder dan 60 ns per patroonpunt . Precompensatie 
van patroon van een gelntegreerde schakeling van circa 
10 10 patroonpunten vergt in dit geval op de huidige 
personal computers derhalve slechts circa 10 minuten. 

15 De uitvinding wordt voorts beschreven in het 

niet-voorgepubliceerde proefschrift met tekst "Proximity 
effects correction in electron beam nanolithography" , 
waarvan de gehele inhoud als hierbij ingelast dient te 
worden beschouwd. 
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CONCLUSIES 

1. Werkwijze voor het bepalen van de voor ver- 
krijging van een gewenst patroon in een deklaag op een 
ondergrond per patroonpositie benodigd geprecompenseerd 
patroon van blootstellingsdoses van een elektronenbundel, 

5 omvattende: 

- het bepalen van de versmeringsf unctie van de 
elektronenbundel ; 

- het met de versmeringsf unctie en het gewenst 
patroon bepalen van het geprecompenseerde patroon, 

10 waarbij het bepalen zodanig wordt uitgevoerd, dat de 

blootstellingsdoses vrijwel uitsluitend positieve waarden 
bevatten en dat de blootstellingsdoses ten opzichte van 
elkaar ten minste in enige mate glad zijn. 

2. Werkwijze volgens conclusie 1, omvattende de 

15 stappen: 

a) het schatten van een regularisatie parame- 
ter; 

b) het met alle patroonpunten van het gewenst 
patroon behalve een bepaald patroonpunt bepalen van een 

2 0 geprecompenseerd patroon; 

c) het met de versmeringsf unctie opnieuw ver- 
smeren van het geprecompenseerd patroon voor het voor- 
spellen van de dosis van het bepaalde patroonpunt; 

d) het voor elk patroonpunt herhalen van stap- 
2 5 pen b en c; 

e) het herhalen van stappen a tot en met d met 
aangepaste regularisatie parameter totdat een eindwaarde 
van een regularisatie parameter bereikt is; 

f) het bepalen van het geprecompenseerde pa- 
30 troon met de eindwaarde van de regularisatie parameter. 

3. Werkwijze volgens conclusie 2, waarbij stap 
b) de volgende iteratieve bepaling omvat : 



14 
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d (j) = d u ' x > + (lTK+kB(D))' 1 K*zU-U r (i) = a-*Td (i) 
met d (0) = 0 en r (0> = a 

waarbij a een vector met als elementen de doses van het 
gewenste patroon, d een vector met de blootstellingsdoses 
5 van het geprecompenseerde patroon, K de versmeringsfunc- 
tie in matrixvorm is, K* de hermitisch geconjugeerde van 
de versmeringsfunctie K is, B een operator, en A. een 
regularisat ieparameter . 

4. Werkwijze volgens conclusie 3, waarbij de 
10 operator B als volgt is gedef inieerd : 

i d tot d tot 

waarin de sommatie plaatsvindt over alle patroonpunten, d. 
het i e element van de vector d is, en d_ de sommatie over 

tot 

15 alle elementen van de vector d voorstelt. 

5. Werkwijze volgens conclusie 2, waarbij in 
stap e) de eindwaarde van de regularisatie parameter de 
regularisatie parameter is waarbij 

20 

waarin N het totaal aantal patroonpunten, a een vector 
met als elementen de doses van het gewenste patroon, d 
een vector met de blootstellingsdoses van het geprecom- 
penseerde patroon, en K de versmeringsfunctie in matrix- 
25 vorm is. 

6. Werkwijze volgens conclusie 2, waarbij in 
stap e) de eindwaarde van de regularisatie parameter de 
regularisatie parameter is waarbij 




30 



T / B E99 / OOO76 



16 



waarin N het totaal aantal patroonpunten, a een vector 
met als elementen de doses van het gewenste patroon, d 
een vector met de blootstellingsdoses van het geprecom- 
penseerde patroon, K de versmeringsf unctie in matrixvorm 
en w^ gedefinieerd is als: 



-1 n 



10 



met a kk de elementen van de matrix A = K (K T K+AL (D) t L (D) ) " 
X K T . en L de Laplace-operator, minimaal is. 

7. Werkwijze volgens een der voorgaande conclu- 
sies, waarbij na stap e) de stappen omvattende van het 
trainen van een neuraal netwerk met behulp van een of 
meer gewenste eerste patronen en de daarbij behorende 
geprecompenseerde patronen. 

8. Werkwijze volgens conclusie 7, waarbij met 
het getrainde neurale netwerk het geprecompenseerde 

15 patroon behorend bij een tweede gewenst patroon te bepa- 
len is. 

9. Werkwijze volgens conclusie 7 en 8, waarbij 
het eerste gewenste patroon een relatief eenvoudig trai- 
ningspatroon is en het tweede gewenste patroon het deel- 

20 patroon van een geintegreerde schakeling is. 

10. Werkwijze volgens conclusie 9, waarbij twee 
of meer deelpatronen tot een samengesteld patroon van de 
geintegreerde schakeling te combineren zijn. 

11. Werkwijze volgens een van de conclusies 7- 
25 10, waarbij het neurale netwerk een radiele basis f unctie 

netwerk is. 

12. Werkwijze volgens een der conclusies 7-11, 
waarbij het neurale netwerk in hardware geimplementeerd 
is . 

30 13. Werkwijze volgens conclusie 12, waarbij het 

neurale netwerk in analoge hardware geimplementeerd is. 



# 

14. Werkwijze volgens een van de voorgaande 
conclusies, waarbij de versmeringsf unctie ten minste is 
opgebouwd uit twee Gaussische functies. 

15. Werkwijze volgens conclusie 14, waarbij aan 
5 de versmeringsf unctie een exponentiele functie is toege- 

voegd . 

16. Werkwijze volgens conclusie 14 of 15, 
waarbij de parameters van de Gaussische functies bepaal- 
baar zijn met behulp van statistische simulaties. 

10 17. Werkwijze volgens conclusie 14 of 15, 

waarbij de parameters van de Gaussische functies bepaal- 

baar zijn door metingen. 

18. Werkwijze voor het bepalen van de voor ver- 

krijging van een gewenst patroon in een deklaag op een 
15 ondergrond per patroonpositie benodigde blootstellingsdo- 

sis van een elektronenbundel , omvattende: 

- het invoeren van een trainingspatroon en de 

daarbij behorende blootstellingsdoses in een neuraal 

netwerk ter training daarvan; 
20 - het bepalen van de weegfactoren van het 

neurale netwerk; 

- het met het neurale netwerk bepalen van de 

voor verkrijging van een gewenst patroon in een deklaag 

op een ondergrond benodigde blootstellingsdoses. 
25 19. Inrichting voor het bepalen van de voor 

verkrijging van een gewenst patroon in een deklaag op een 

ondergrond per patroonpositie benodigde blootstellingsdo- 

sis van een elektronenbundel, omvattende elektronische 

schakelmiddelen voor het implementeren van een neuraal 
3 0 netwerk met weegfactoren die bepaald zijn volgens een der 

voorgaande conclusies . 

20. Geintegreerde schakelingen die zijn ver- 

vaardigd met de inrichting van conclusie 19 of volgens de 

werkwijze van een der conclusies 1-18. 
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UITTREKSEL 

De onderhavige uitvinding heeft betrekking op 
een werkwijze voor het bepalen van de voor verkrijging 
van een gewenst patroon in een deklaag op een ondergrond 
per patroonpositie benodigd geprecompenseerd patroon van 
5 blootstellingsdoses van een elektronenbundel, omvattende: 

- het bepalen van de versmeringsf unctie van de 
elektronenbundel ; 

- het met de versmeringsf unctie en het gewenst 
patroon bepalen van het geprecompenseerde patroon, 

10 waarbij het bepalen zodanig wordt uitgevoerd, dat de 

blootstellingsdoses vrijwel uitsluitend positieve waarden 
bevatten en dat de blootstellingsdoses ten opzichte van 
elkaar ten minste in enige mate glad zijn. 
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CLAIMS 

1. Method for determining the precompensated 
pattern of exposure doses of an electron beam required 
per pattern position to obtain a desired pattern in a 
coating on a substrate, comprising of: 

5 - determining the smearing function of the 

electron beam; 

- determining the precompensated pattern with - 
the smearing function and the desired pattern, wherein 
the determination is performed such that the exposure 
10 doses contain almost exclusively positive values and that 
the exposure doses are at least to some extent smooth 
relative to each other. 

2. Method as claimed in claim 1, comprising the 

steps of: 

15 a) estimating a regularization parameter; 

b) determining a precompensated pattern with 
all pattern points of the desired pattern with the 
exception of a determined pattern point; 

c) smearing the precompensated pattern again 

2 0 with the smearing function in order to predict the dose 

of the determined pattern point; 

d) repeating steps b and c for each pattern 

point; 

e) repeating steps a to d with adapted 
25 regularization parameter until a final value of a 

regularization parameter is obtained; 

f ) determining the precompensated pattern with 
the final value of the regularization parameter. 

3. Method as claimed in claim 2, wherein step 

3 0 b) comprises the following iterative definition: 

d (D = d^ + UCK+XBiD) ) -i^r (i - 1J r (i) = a-Kd U) 



with d t0) = 0 and r< 0) = a 
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wherein a is a vector with the doses of the desired 
pattern as elements, d is a vector with the exposure 
doses of the precompensated pattern, K is the smearing 
function in matrix form, K" is the Hermitian conjugate of 
5 the smearing function K, B is an operator and A a 
regularization parameter. 

4. Method as claimed in claim 3, wherein the 
operator B is defined as follows : 

B{D)-Y, <-^->ln(-^-) 



10 



d toc d coc 



in which the summation takes place over all pattern 
points, d. is the i ch element of the vector d, and d 

cot 

represents the summation over all elements of the vector 
d. 

15 5 - Method as claimed in claim 2, wherein the 

final value of the regularization parameter in step e) is 
the regularization parameter wherein 



jff: (a* - Ud k (k)l k ) 



wherein N is the total number of pattern points, a is a 
20 vector with the doses of the desired pattern as elements, 

d is a vector with the exposure doses of the 

precompensated pattern and K the smearing function in 

matrix form. 

6. Method as claimed in claim 2, wherein the 
25 final value of the regularization parameter in step e) is 

the minimal regularization parameter wherein 



i N 



a k - [Kd k (X) } k ) 2 w kk {\) 
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wherein N is the total number of pattern points, a is a 
vector with the doses of the desired pattern as elements, 
d is a vector with the exposure doses of the 
precompensated pattern, K is the smearing function in 
5 matrix form and w klt is defined as: 




with a kk the elements of the matrix A = K (K T K+AL (D) T L (D) ) *- 
1 K T and L the Laplace operator. 

7. Method as claimed in any of the foregoing 
claims, wherein after step e) the step is performed of 

10 training a neural network using one or more desired first 
patterns and the associated precompensated patterns. 

8. Method as claimed in claim 7, wherein the 
precompensated pattern associated with a second desired 
pattern can be determined with the trained neural 

15 network. 

9 . Method as claimed in claims 7 and 8 , wherein 
the first desired pattern is a relatively simple training 
pattern and the second desired pattern is the partial 
pattern of an integrated circuit. 

20 10. Method as claimed in claim 9, wherein two 

or more partial patterns can be combined into a composite 
pattern of the integrated circuit. 

11. Method as claimed in any of the claims 7- 

10, wherein the neural network is a radial base function 
25 network-. 

12. Method as claimed in any of the claims 7- 

11, wherein the neural network is implemented in 
hardware . 

13. Method as claimed in claim 12, wherein the 
3 0 neural network is implemented in analog hardware. 

14 . Method as claimed in any of the foregoing 
claims, wherein the smearing function is made up of at 
least two Gaussian functions . 
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15. Method as claimed in claim 14, wherein an 
exponential function is added to the smearing function. 

16. Method as claimed in claim 14 or 15, 
wherein the parameters of the Gaussian functions can be 

5 determined using statistical simulations. 

17. Method as claimed in claim 14 or 15, 
wherein the parameters of the Gaussian functions can be 
determined by measurements . 

18. Method for determining the exposure dose of 
10 an electron beam required per pattern position to obtain 

a desired pattern in a coating on a substrate, comprising 
of: 

- inputting a training pattern and the 
associated exposure doses into a neural network for 

15 training thereof; 

- determining the weighting factors of the 
neural network; 

- determining with the neural network the 
exposure doses required to obtain a desired pattern in a 

2 0 coating on a substrate. 

19. Device for determining the exposure dose of 
an electron beam required per pattern position to obtain 
a desired pattern in a coating on a substrate, comprising 
electronic circuit means for implementing a neural 

25 network with weighting factors determined as claimed in- 
any of the foregoing claims. 

20. Integrated circuits manufactured with the 
device of claim 19 or according to the method of any of 
the claims 1-18. 
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METHOD AND DEVICE FOR CORRECTING PROXIMITY EFFECTS 

The present invention relates to a method and 
device for determining the exposure dose of an electron 
beam required per pattern position to obtain a desired 
pattern in a coating. 
5 In the manufacture of the latest generations of 

integrated circuits use is preferably made of focused 
electron beams in lithographic processes instead of 
making use of the usual optical lithographic techniques, 
since these latter techniques are subject to limitations 

\0 in terms of resolution as a result of diffraction of the 
used laser light. The resolution of the integrated cir- 
cuit obtained with such electron beam lithography is 
greater, although it is limited by scattering of the 
electrons in the coating. Methods are known for minimi - 

15 zing scattering effects or compensating therefor in 
advance and thereby improving the resolution of the 
obtained integrated circuits. 

The known methods have the drawback however 
that scattering effects can themselves only be minimized 

2 0 to a limited degree, while advance compensation according 

to the known method requires many calculations and there- 
fore needs a long calculation time. For the manufacture 
of integrated circuits for instance a very large number 
Of pattern points, often in the order of magnitude of 10 10 
25 pattern points, must be "written" , while the number of 
calculations required for this purpose amounts to a 
multiple thereof. As a result a practically real time 
precompensation for the smearing or blurring effects 
cannot be implemented. 

3 0 The object of the present invention is to 

obviate this drawback and also provide additional advan- 
tages . 

The present invention therefore relates to a 
method for determining the precompensated pattern of 
3 5 exposure doses of an electron beam required per pattern 
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position to obtain a desired pattern in a coating on a 
substrate, comprising of: 

- determining the smearing function (blur 
function) of the electron beam; 
5 - determining the precompensated pattern with 

the smearing function and the desired pattern, wherein 
the determination is performed such that the exposure 
doses contain almost exclusively positive values and that 
the exposure doses are at least to some extent smooth 

10 relative to each other. 

Since a negative value for the exposure doses 
of an electron beam has no physical significance and 
cannot therefore be realized, the determination of the 
exposure doses of the precompensated pattern is performed 

15 such that it assumes (almost) exclusively positive valu- 
es. A smooth solution is furthermore obtained since 
strong oscillations in the smearing function have no 
physical basis but are caused solely by mathematical 
instability of the calculations. 

2 0 In a preferred embodiment of the invention the 

method comprises the steps of : 

a) estimating a regularization parameter; 

b) determining a precompensated pattern with 
all pattern points of the desired pattern with the excep- 

2 5 tion of a determined pattern point; 

c) smearing the precompensated pattern again 
with the smearing function in order to predict the dose 
of the determined pattern point; 

d) repeating steps b and c for each pattern 

3 0 point; 

e) repeating steps a to d with adapted regula- 
rization parameter until a final value of a regularizati- 
on parameter is obtained; 

f) determining the precompensated pattern with 
35 the final value of the regularization parameter. 

According to a further embodiment of the inven- 
tion step b) comprises the following iterative definiti- 
on : 
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with d (0) = 0 and r (0) = a 

in which a is a vector with the doses of the desired 
pattern as elements, d is a vector with the exposure 
5 doses of the precompensated pattern, K is the smearing 
function in matrix form, K* is the Hermitian conjugate of 
the smearing function K, B is an operator and X a 
regularization parameter. 

According to a further embodiment of the 
10 invention operator B is defined as follows: 

B(D) ( 7 ^-)in(-A-) 

i "tot "tot 

in which the summation takes place over all pattern 
points, d i is the i th element of the vector d, and d 
represents the summation over all elements of the vector 
d. 

According to a further embodiment of the 
invention the final value of the regularization parameter 
in the above mentioned step e) is the regularization 
parameter wherein 

jri: <«* - iKd k a)] k )z 

in which N is the total number of pattern points, a is a 
vector with the doses of the desired pattern as elements, 
d is a vector with the exposure doses of the 

2 5 precompensated pattern and K the smearing function in 

matrix form. 

According to a further embodiment of the 
invention the final value of the regularization parameter 
in step e) is the minimal regularization parameter 

3 0 wherein 
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±£ (a k - [Kd*(\)] k )*w kk a) 

in which N is the total number of pattern points, a is a 
vector with the doses of the desired pattern as elements, 
d is a vector with the exposure doses of the 
precompensated pattern, K is the smearing function in 
5 matrix form and w.. is defined as: 




with a kk the elements of the matrix A = K (K T K+A.L (D) T L (D) ) " 
L K T and L the discrete Laplace transformation. 

According to a further embodiment of the 
invention after step e) the step is performed of training 

10 a neural network using one or more desired first patterns 
and the associated precompensated patterns . 

According to a further embodiment of the 
invention the precompensated pattern associated with a 
second desired pattern can be determined with the trained 

15 neural network, wherein in a further embodiment the first 
desired pattern is a relatively simple training pattern 
and the second desired pattern is the partial pattern of 
an integrated circuit, and wherein in a further 
embodiment two or more partial patterns can be combined 

20 into a composite pattern of the integrated circuit. 

By determining the associated precompensated 
pattern of exposure doses for a known desired pattern, 
which is preferably simple, and then establishing the 
relation between the weighting factors of a neural 

25 network, is ensured that for a second desired pattern, 
which may be complicated, obtaining the relation between 
this pattern and the associated exposure doses is 
determined in very efficient and rapid manner by the 
neural network. The first pattern is generally a 

30 relatively simple training pattern, while the second 
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pattern is for instance the pattern of a very complicated 
integrated circuit . 

In a preferred embodiment of the invention the 
above stated neural network is implemented in hardware, 
5 whereby determining of the relation between a pattern and 
the exposure dose associated therewith is performed in 
more rapid manner, for instance within 6 0 ns for each 
pattern point and within 10 minutes for a pattern of 10 10 . 
According to a preferred embodiment of the 

10 invention the smearing function is made up of at least 
two Gaussian functions, to which an exponential function 
is optionally added. Parameters of the Gaussian functions 
and optionally the exponential function can be determined 
by means of statistical simulation of the system of 

15 electron beam transmitting equipment and the relevant 
coating and substrate of the integrated circuit for 
manufacture . 

In another embodiment of the invention 
parameters are determined by performing measurements on 

20 the system of electron beam transmitting equipment and 
the relevant coating with substrate. 

The present invention also relates to a device 
for determining the exposure dose of an electron beam 
required per pattern position to obtain a desired pattern 

25 in a coating on a substrate, comprising an electronic 
circuit for implementing said neural network with 
weighting factors determined in the above stated manner. 

The invention will be elucidated hereinbelow 
with reference to a preferred embodiment thereof, wherein 

3 0 use will be made of the annexed drawings, in which: 

- figure 1 shows a schematic overview of a 
preferred embodiment of a device according to the 
invention; 

- figures 2a- 2c show a schematic overview of 
35 the determination of a precompensated pattern of 3x3 

pattern points; 

- figure 3 shows a desired training pattern of 
256x256 pattern points; 
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ft - figure 4 shows the training pattern of figure 
3 after smearing; 

- figure 5 shows a graph in which for the 
training pattern of figure 3 the prediction error is 

5 plotted as a function of the chosen regularization 
parameter; 

- figure 6 shows the training pattern of figure 
3 after precompensation; 

- figure 7 shows the precompensated pattern of 
10 figure 6 after smearing; and 

- figure 8 is a schematic representation of a 
neural network for determining precompensated patterns. 

In an arrangement of equipment for transmitting 
an electron beam and a substrate 1 with coating 2 for 

15 processing, a beam of electrons 3 is directed at a 

position or pattern point of a coating 2 on a substrate 
1 . The interaction of the incident electron beam 3 with 
the coating or resist film 2 and the underlayer or 
substrate 1 results in a scattering of the electrons in 

20 coating 2 which causes smearing or proximity effects. 

When for instance a primary electron penetrates into the 
coating, a part of its energy is transferred to electrons 
of the atoms of the coating, which causes ionization or 
excitation thereof. A collision between electrons with a 

25 large transfer of energy generates a secondary electron 
which generally has a direction of movement perpendicular 
to that of a primary electron. 

Smearing effects in electron beam lithography 
relate more generally to the process whereby the 

3 0 resolution of the exposed pattern is reduced by primary 
electron scattering (forward scattering) and secondary 
electron excitation (backward scattering) in the coating 
and the substrate of an integrated circuit for 
manufacture. Sharp features such as angles in the desired 

35 pattern are rounded, line thicknesses and interspaces are 
modified and in particular extreme cases some features 
even disappear completely or are merged in incorrect 
manner with adjacent features. 
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The smearing effects or proximity effects can 
be described by a smearing function, which shows the 
relation between on the one hand the exposure doses of a 
determined pattern point of a pattern for manufacturing 
5 in the coating and on the other the doses actually 
absorbed by this pattern point and adjacent pattern 
points. The effect of the smearing is thus established in 
the smearing function. 

Assuming that exposure and smearing are 

10 linearly and spatially invariant and that for a numeric 
solution a discrete representation is preferred, the 
above can be expressed in matrix form as follows: A = KD, 
in which A is a column vector of which each element a is 
the total energy dose which is actually absorbed in the 

15 associated pattern point, K is a smearing matrix of which 
each mn th element is the portion of the energy dose which 
is absorbed in pattern point m from a unit -exposure dose 
supplied to pattern point n, and D is a column vector 
made up of elements d. which represent the exposure doses 

20 generated per pattern point by the electron beam 

equipment. Since the smearing effect is unavoidable, it 
is best to adapt the exposure doses d. of the different 
pattern points such that the dose a i actually absorbed in 
a pattern point is such that the desired pattern is still 

25 obtained. 

This so-called precompensat ion of the exposure 
dose of the electron beam can be performed according to 
the prior art by determining the inverse of the smearing 
matrix K. There are many ways of generally inverting a 

3 0 matrix. However, these methods often take no account of 
physical limitations, such as in this case for instance 
those of the electron beam transmitters. No negative 
exposure doses for instance are thus possible. A further 
drawback of such inversion methods is that the inverted 

35 matrix has many oscillations. In addition, for inversion 
of the smearing matrix for a partial pattern of for 
instance 256x256 pattern points the inversion of a 
smearing matrix with dimensions of 65536x65536 has to be 
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calculated, which requires an enormous amount of 
calculating time . 

Figures 2a to 2c show a desired pattern (A) . 
The pattern is built up of 9 pattern points a. wherein i 
5 varies from 1 to 9 . This desired pattern must be 

precompensated in order to be able to provide the desired 
pattern after exposure to the smearing electron beam, 
i.e. the values of d., with i varying from 1 to 9 , have to 
be determined . 
10 The precompensated pattern is first of all 

determined making use of the doses a with i from 2 to 9, 
wherein pattern point 1 is not therefore taken into 
account (figure 2a) . This precompensated pattern is 
determined on the basis of the following expression: 

d d) = d (1 - 1) + (K*K+XB(D))- 1 K*r il - 1) z U) = a-Kd U) 

15 with d (0) = 0 and r< 0) = a < 

wherein a is a vector with the doses of the desired 
pattern as elements, d is a vector with the exposure 
doses of the precompensated pattern , K is the smearing 
function in matrix form, K* is the Hermitian conjugate of 

20 smearing function K, B is an operator and X a 

regularization parameter. The value of the regularization 
parameter can be chosen at random, in this case for 
instance X = 0 . 

The operator B imposes a limitation and can be 

25 defined as follows: 

B(D) ( J^)i n( .A-) 

in which the summation takes place over all pattern 
points, d is the i th element of the vector d, and d 

1 tot 

30 represents the summation over all elements of the vector 
d. 

The thus determined precompensated pattern is 
then smeared once again on the basis of the known 
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smearing function, whereby the predicted dose Kd of 
pattern point 1 is determined. 

The above procedure is then repeated 
successively (figures 2b and 2c) for the second to ninth 
5 pattern point (i=2,..,9), wherein all pattern points with 
the exception of one pattern point are used each time. 

On the basis of the above results, the least 
squares prediction error over all pattern points is 
determined, which will be further explained later. 

10 The above procedure is subsequently repeated 

with different values for the regularization parameter X , 
for instance X 2 = 0.001, X= 0.002 etc. The regularization 
parameter is eventually chosen wherein the least squares 
prediction error over all pattern points is minimal. This 

15 regularization parameter is then chosen as the optimal 
regularization parameter A opc . The final precompensated 
pattern is then determined on the basis of this optimal 
regularization parameter A . 

For this purpose the minimum is determined of 

2 0 the expression: 

j-j; <«* - iKd*(X)] k )>w kk a) 

1V Jc=l 

in which N is the total number of pattern points, a is a 
vector with the doses of the desired pattern as elements, 
25 d is a vector with the exposure doses of the 

precompensated pattern, K is the smearing function in 
matrix form and w^ is defined as: 

iv J=l 

with a kk the elements of the matrix A = K (K T K+AL (D) t L (D) ) " 
1 K T and L the Laplace operator. 
30 The smearing function resulting from forward 

scattering and backward scattering of the electrons of 
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the electron beam can be determined in different ways. It 
can be determined on the basis of measurements of the 
impulse" response of the equipment for transmitting the 
electron beam on a test object. The smearing function can 
5 also be determined using diverse Monte Carlo techniques. 
In the first method of determination all physical aspects 
of the equipment used are taken into account. In the 
latter mentioned method of determination only a model of 
the reality is used, although the determination is 
10 however easier to perform without requiring extensive 
measurements . 

Gaussian functions are preferably used as 
approximation for the smearing functions determined in 
any of the above described methods. The smearing function 
15 is in this case "fitted" for instance with a scattering 
fit model of a double Gaussian function (for both forward 
and backward scattering properties of the electrons) , a 
triple Gaussian function or a double Gaussian function 
with a decreasing exponential function. The choice of the 

2 0 scattering fit model depends on the dimensions of the 

components to be distinguished in the test object 
(resolution) . At dimensions smaller than 100 nm the 
choice hereof becomes critical: at such small dimensions 
the triple Gaussian functions or double Gaussian 
25 functions with decreasing exponential function are 

recommended. A smearing function with double Gaussian 
function can be described with 3 parameters, while the 
other two stated scattering fit models can be described 
with 4 parameters, which implies a great reduction in the 

3 0 quantity of data for processing. 

Figure 3 shows a desired pattern of 256x256 
pattern points. Smearing with a smearing function in the 
form of a double Gaussian function with a= 50 nm, S= 3.45 
and n = 1.36 produces the smeared pattern of figure 4. It 
35 is clearly visible that much detail in the pattern has 
been lost, which means a limitation in the resolution to 
be obtained of the pattern for manufacture. Application 
of the method according to the invention produces an 
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optimal regularization parameter of A. = 0.07042, which is 
shown in* figure 5, in which the error in the pattern is 
minimal at this value of X. The precompensated pattern 
calculated with this value of k and the associated 
5 smeared pattern are shown respectively in figures 6 and 
7. Comparison of the results of figure 7 with those of 
figure 3 shows that the precompensation of the pattern 
with a desired pattern produces a smeared pattern with a 
greatly improved resolution. Components for 

10 distinguishing with dimensions of less than 100 nm, for 
instance in an integrated circuit, can hereby be 
realized. A comparison of the results of the method 
described herein with those of other correction methods 
is shown in table 1. The degree of error of the 

15 correction methods is defined here as the summation of 
the difference between the calculated exposure doses and 
the ideal precompensated exposure doses divided by the 
number of pattern points. 



Correction method 


degree of error in % 


Uncorrected 


10 .2' % 


Truncating 


10 .2 % 


Shifting and scaling 


12.2 % 


- Present method 


4.9% 



25 

From the above can be seen that the present method of 
determining a precompensated pattern produces by far the 
smallest degree of error compared to the other usual 
methods . 

3 0 The precompensated pattern and the desired 

pattern are subsequently used as training set or training 
patterns for a neural network. A part of such a network 
is shown schematically in figure 8 en is represented by 
the expression 
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. i.e. the dose a i is expressed in a set of 9 basic 
functions h. , in this case radial functions. 

After training of the neural network a 
5 precompensated pattern can be determined for another 
random desired pattern in very rapid manner. A random 
pattern can for instance be a pattern of 512 by 512 
pattern points forming a partial pattern of an integrated 
circuit. Various partial patterns can then be combined 
10 (clustered) to form one pattern which comprises the whole 
integrated circuit or at least a part thereof. 

The above described neural network can be 
implemented in hardware, and preferably in analog 
hardware since the calculating speed of neural networks 
15 implemented in this manner is very great. The calculating 
time for precompensation of a pattern thus amounts to 
less than 60ns per pattern point. Precompensation of a 
pattern of an integrated circuit of about 10 10 pattern 
points requires in this case only -about 10 minutes on 
2 0 present personal computers. 

The invention is further described in the non- 
prepublished doctoral thesis with the title "Proximity 
effects correction in electron beam nanolithography " , the 
entire content of which should be deemed as interpolated 
25 herein. 
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CLAIMS 

1. Method for determining the precompensated 
pattern of exposure doses of an electron beam required 
per pattern position to obtain a desired pattern in a 
coating on a substrate, comprising of: 

5 - determining the smearing function of the 

electron beam; 

- determining the precompensated pattern with 
the smearing function and the desired pattern, wherein 
the determination is performed such that the exposure 
10 doses contain almost exclusively positive values and that 
the exposure doses are at least to some extent smooth 
relative to each other. 

2. Method as claimed in claim 1, comprising the 

steps of: 

15 a) estimating a regularization parameter; 

b) determining a precompensated pattern with 
all pattern points of the desired pattern with the 
exception of a determined pattern point; 

c) smearing the precompensated pattern again 
20 with the smearing function in order to predict the dose 

of the determined pattern point; 

d) repeating steps b and c for each pattern 

point ; 

e) repeating steps a to d with adapted 
25 regularization parameter until a final value of a 

regularization parameter is obtained; 

f) determining the precompensated pattern with 
the final value of the regularization parameter. 

3. Method as claimed in claim 2, wherein step 
30 b) comprises the following iterative definition: 

d (j) = d^ + iK'K+XBiD) ) -ijrr^- 1 * r (J) = a-Kd i2) 

with d (0) = 0 and r< 0) = a 
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wherein a is a vector with the doses of the desired 
pattern as elements, d is a vector with the exposure 
doses of the precompensated pattern, K is the smearing 
function in matrix form, K* is the Hermitian conjugate of 
5 the smearing function K, B is an operator and X a 
regularization parameter. 

4 . Method as claimed in claim 3, wherein the 
operator B is defined as follows: 

B(D) -£ (-A-)in(-^) 

i a tot a tot 

10 

in which the summation takes place over all pattern 
points, d i is the i th element of the vector d, and d 
represents the summation over all elements of the vector 
d. 

15 5. Method as claimed in claim 2, wherein the 

final value of the regularization parameter in step e) is 
the regularization parameter wherein 

lv ic=l 

wherein N is the total number of pattern points, a is a 
20 vector with the doses of the desired pattern as elements, 

d is a vector with the exposure doses of the 

precompensated pattern and K the smearing function in 

matrix form. 

6. Method as claimed in claim 2, wherein the 
25 final value of the regularization parameter in step e) is 

the minimal regularization parameter wherein 
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wherein N is the total number of pattern points, a is a 
vector with the doses of the desired pattern as elements, 
d is a vector with the exposure doses of the 
precompensated pattern, K is the smearing function in 
5 matrix form and w kk is defined as: 

iV j=i 

with a kJc the elements of the matrix A = K (K T K+AL (D) t L (D) ) * 
1 K T and L the Laplace operator. 

7 . Method as claimed in any of the foregoing 
claims, wherein after step e) the step is performed of 

10 training a neural network using one or more desired first 
patterns and the associated precompensated patterns. 

8. Method as claimed in claim 7, wherein the 
precompensated pattern associated with a second desired 
pattern can be determined with the trained neural 

15 network. 

9. Method as claimed in claims 7 and 8, wherein 
the first desired pattern is a relatively simple training 
pattern and the second desired pattern is the partial 
pattern of an integrated circuit. 

20 10. Method as claimed in claim 9, wherein two 

or more partial patterns can be combined into a composite 
pattern of the integrated circuit. 

11. Method as claimed in any of the claims 7- 

10, wherein the neural network is a radial base function 
25 network. 

12. Method as claimed in any of the claims 7- 

11, wherein the neural network is implemented in 
hardware . 

13. Method as claimed in claim 12, wherein the 
3 0 neural network is implemented in analog hardware. 

14 . Method as claimed in any of the foregoing 
claims, wherein the smearing function is made up of at 
least two Gaussian functions. 
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15. Method as claimed in claim 14, wherein an 
exponential function is added to the smearing function. 

16. Method as claimed in claim 14 or 15, 
wherein the parameters of the Gaussian functions can be 

5 determined using statistical simulations. 

17. Method as claimed in claim 14 or 15, 
wherein the parameters of the Gaussian functions can be 
determined by measurements. 

18. Method for determining the exposure dose of 
10 an -electron beam required per pattern position to obtain 

a desired pattern in a coating on a substrate, comprising 
of: 

- inputting a "training pattern and the 
associated exposure doses into a neural network for 

15 training thereof; 

- determining the weighting factors of the 
neural network ; 

- determining with the neural network the 
exposure doses required to obtain a desired pattern in a 

2 0 coating on a substrate. 

19. Device for determining the exposure dose of 
an electron beam required per pattern position to obtain 
a desired pattern in a coating on a substrate, comprising 
electronic circuit means for implementing a neural 

25 network with weighting factors determined as claimed in 
any of the foregoing claims . 

20. Integrated circuits manufactured with the 
device of claim 19 or according to the method of any of 
the claims 1-18. 
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